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ABSTRACT 

Context. Within the collision growth scenario for planetesimal formation, the growth step from centimetre-sized pre-planetesimals 
to kilometre-sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be 
halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount 
of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of 
pre-planetesimal collisions is required and has to be presented in a suitable and precise format. 

Aims. We wish to develop a new classification scheme broad enough to encompass all events with sticking, bouncing, and fragmen- 
tation contributions, accurate enough to capture the important collision outcome nuances, and at the same time simple enough to be 
implementable in global dust coagulation simulations. We furthermore wish to demonstrate the reliability of our numerical smoothed 
particle hydrodynamics (SPH) model and the applicability of our new collision outcome classification to previous results as well as 
our simulation results. 

Methods. We propose and apply a scheme based on the quantitative aspects of four fragment populations: the largest and second 
largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we 
adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. On the basis of laboratory benchmark ex- 
periments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation 
behaviour of macroscopic highly porous Si02 dust aggregates. 

Results. We demonstrate that previous attempts to map collision data were much too oriented on qualitatively categorising into stick- 
ing, bouncing, and fragmentation events. Intermediate categories are found in our simulations that are difficult to map to existing 
qualitative categorisations. We show that the four-population model encompasses all previous categorisations and in addition allows 
for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, 
and filling factor. In addition, the numerical porosity model successfully passes another benchmark test: the correct simulation of the 
entire list of collision outcome types yielded by laboratory experiments. As a demonstration of the applicability and the power of 
the four-population model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of 
intermediate porosity aggregates. 

Key words, accretion, accretion discs - hydrodynamics - methods: data analysis - methods: numerical - planets and satellites: forma- 
tion - protoplanetary discs 

1. Introduction tion, drift, and bouncing barrier (e.g. IZsom et aP l2010l) . The 

most serious barrier is the fragmentation barrier: with increas- 
The primary way in which rocky planets are thought to form size, the relative velocities between pre-planetesimals in- 

is by core accretion in protoplanetary discs consisting of gas crease and potentially lead to catastrophic disruption in mu- 
and dust. Interaction between these components induces rela- tual colUsions. Often a velocity threshold of 1 ms ' for disrup- 
H tive velocities and consequential mutu al collisions among the tiye events is assumed (iBlum & MunchllT993l) . The threshold is 
9 initially micron-sized dust grain s (e.g. |Weidenschining||l977bi mosdy assumed to be as independent of other par ameters such 
IWeidenschiUing & Cuzzi| \l993- In this way, millimetre- to as_mrosity, mass ratio, and impact parameter (e.g. |GuttleretalJ 
centimetre-sized fluffy dust aggregates grow by a simple hit- ^^mh . With the assumption of thi s fragmentation thre shold and 
and-stick mechanism (e.g. Blum & Wurm 20Q&). The growth power-law distributed fragments, iBrauer et al.l (l2008ll find that 
mechanism in the pre-planetesimal regime (centimetre to kilo- dust coagulation is halt ed at centimetre or even milhmetre sizes, 
metre) is unclear and is addressed in this paper Once a suflicient iTeiser & Wuriiil (120091) proposed a growth model that is based 
population of kilomene-sized objects has formed, gravitation- mainly on the accretion of fragments < 1 mm by larger objects 
assisted accretion ensures the final assembly of full-sized plan- and therefore requires some amount of disruptive coUisions. 
ets (e.g. lGoldreich et al...2004() . The planetary building blocks in 

this gravity-dominated regime are referred to as planetesimals. The second obstacle to planetesii nal formation is flie ra - 

The growth step from pre-planetesimals to planetesimals is dial drift of larger pre-planetesimals (IWeidenschillind [l 977ah . 
not well-consnained and the subject of extensive numerical and Owing to its pressure support, the gas in the disc rotates at a sub- 
experimental effort. The difficulties to overcome in this step Keplerian velocity, whereas the solid material lacks this pressure 
can be characterised with three types of barriers: fragmenta- support and tends to rotate at a Keplerian velocity. As a con- 
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sequence, the solid objects feel a headwind that causes an in- 
ward drift. For the minimum-mass solar nebula, it takes roughly 
a century for metre sized objects to drift from 1 AU into the 
star. In contrast, for centimetre- and kilometre-sized objects it 
takes ~ 10^ yr. This is because small objects are strongly cou- 
pled to the motion of the gas, whereas large objects are almost 
decoupled. As a consequence, in the intermediate regime, metre- 
sized objects are quickly lost by accretion onto the host star or 
by photoevaporation in the hot zones close to it. Collective ef- 
fects in the midplane of the protoplanetary disc might dimin- 
ish the problem^ (e.g ICuzzi et al.lil993L iDominik et al] l2007l 
IWeidenschilling||20iq) . 

On the basis of empirical dat a showing the exis- 
tence of rebound in dust collisions (iBlum & Munchl I1993L 



HeiB elmann et alj|2007l 'Langko wski et al.H2008L IWeidling et all 
2009 . Guttler et al.l2010 ). Zsom et al.l(l2010h ^oduced a possi- 
ble new obstacle to growth: the bouncing barrier, where collision 
growth is halted at centimetre-sized pebbles. In the collisions, 
aggregates only become compacted in mutual collisions but do 
not grow any further However, recent studies (Geretshauser et 
al., in prep., Wada et al., in prep.) indicate that the collision pa- 
rameter space go verned by bouncing may have been estimated 
to be too large bv lGiittler et al.l (1201 Ol) and that rebound is not as 
frequent as they inferred. 

However, frag mentation is also necessary in discs. 
Observations (e.g. iNatta et al.ll2b07h suggest that sub-millimetre 
sized dust and millimetre- to centimetre-sized pebbles are 
present for as long as about 1 0^ yr in T Tauri discs. However, 
iDullemond & DominikI (l2005h find that without fragmentation 
the disc is depleted of small grains within lO^* yr. This is 
commonly referred to a s the grain re t ention proble m. With 
coagulation simulations, iBimstiel et al.l (l2009l l20I0h tried to 
match the fragment distribution of disruptive events with the 
observational data. 

To summarise, the formation of planetesimals requires 
the right amount of sticking, bouncing, and fragmentation to 
be consistent with observations. Therefore, collisions of pre- 
planetesimals have to be investigated as thoroughly as possible 
and their outcome has to be mapped as precisely as possible tak- 
ing into account all the relevant parameters such as initial poros- 
ity, collision partner size, impact vel ocity, mass ra t io, an d rota- 
tion. To classify collision outcomes, iGiittler et al.l (1201 Ol) com- 
piled 19 experiments and mapped them to four types of stick- 
ing, two types of bouncing, and three types of fragmentation. 
However, owing to experimental restrictions they covered only 
small parts of the relevant parameter space. Many of their find- 
ings were not deduced from collisions between porous dust ag- 
gregates, but from dust collisions with a solid object. Collisions 
between aggregates larger than the decimetre size were not pos- 
sible because of the limitations of their experimental apparatus. 
In addition, not all experiments could be carried out in protoplan- 
etary disc conditions, i.e., in a vacuum and under microgravity. 

An alternative to experiments is to study pre-planetesimal 
collisions numerically. To explore the region where experimen- 
tal data is lacking, we perform simulations of pre-planetesimal 
collisions with our solid-body smoothed particle hydrodynam- 
ics (S PH) code, which was previously used by Schiifer et al] 
(l2007h to emphasise the importance of precise material param- 
eters a nd a thorough calibration of the implemented porosity 
model. iGeretshauser et al] (l20I0h expanded the porosity model 
and both calibrated and tested it extensively to simulate pre- 
planeta ry dust material ba sed on laboratory benchmark exper- 
iments dGiittler et al.ll2009l) . The code and calibration procedure 
is briefly described in Sec.|2] In Sec.|3] we show that the code can 



furthermore reprod uce all sticking, boun cing, and fragmentation 
types proposed bv IGiittler et al.l (I2OIOI) . However, we find that 
their categorisation can also introduce unnecessary complexity 
and on some occasions may lack the required accuracy. It is pos- 
sible that in a collision more than one process of sticking, bounc- 
ing, and fragmentation can take place and qualitative models do 
not make it clear whether the overall growth is positive, negative, 
or neutral. Therefore, we propose a new, simpler but at the same 
time more quantitative model for mapping collision data that 
is presented in Sec. |4] To show the applicability of this model, 
we present the first results of simulations of collisions between 
macroscopic objects consisting of realistic pre-planetesimal ma- 
terial in Sec. |5] We finally summarise our findings and discuss 
future work in Sec. |6] 



2. Code and calibration 



For the simulations pres ented in this paper, we utilise the paral- 
lel SPH code paras ph dHipp & Rosenstielll2004l ISchafeil^OOl 
Is chafer etal.1 [20071). which was calibrat ed to simulate porous 
dust material ( Geretshauser et akl I2OIOI) . This section briefly 
summarises the numerical model and calibration procedure. 



2. 1 . SPH and porosity model 

For collision simulations, SPH as a meshless Lagrangian par- 
ticle method is superior to grid-based methods. This is be- 
cause the continuous fluid or solid is divided into inter- 
acting mass packages, which serve as numerical sampling 
points and form a natural frame of referen ce for deforma- 
tion and fragmentati on (see reviews, e.g. by lMonaghanll2005l 
and iRosswo j l2009t) . Originally developed for astrophysical 
fluid flows , SPH has also been extended to model solid ma- 
terial (e.g. 'Liberskv & Petschek' '1991', 'Benz & Asphauglll994i 
iRandles & Liberskv„1996. Liberskv et al.,i,1997i) . 

Within the SPH scheme, the equations of continuum me- 
chanics are solved in their Lagrangian form. In this framework, 
the continuity and momentum equations are given by 



(1) 
(2) 



where the Greek indices indicate the spatial components and the 
Einstein sum convention is applied. The quantities p and v are 
the density and velocity, respectively, and ctq,^ is the stress tensor 
which is defined as 
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where p is the pressure, which accounts for pure hydrostatic 
compression or tension, and Sap is the deviatoric stress tensor, 
which represents pure shear. For the time evolution of the devia- 
toric stress tensor in the elastic regime, we adopt the approach of 
iBenz&Asphaud ( 119941) . which uses the Jaumann rate form and 
follows Hooke's law. The elastic hydrostatic pressure is mod- 
elled as part of the porous equation of state (Eq. |8}. 

To model sub-resolution porosity, we relate the filling factor 
(f) to the continuous quantities according to 



Ps 



(4) 
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where ps = 2000 kg m ^ is the dens ity of the matrix material 
of Si02 (e.g. lBlum & S chi-at3ler"2004. Table [T]shows the mate- 
rial parameters of t he cal ibrated porosity model determined by 
iGeretshauser et all (l2010h. To model th e plastic behaviour, we 
modified the approach bv lSironol (l2004h . who models plasticity 
with filling-factor-dependent strength quantities. The compres- 
sive strength 2(0), the tensile strength and shear strength 
Y{(p) represent transition thresholds between the ela stic and plas- 
tic reg im e. In particular, we adopt t he relations (see lGiittler et al.l 
120091 and lGeretshauser et"ani2010l) 



2(0) = p„ 



fore 



A In 10 



(5) 



bmin+s < (f) < 0max and s = 0.005. The quantities 0n,in = 0.12 
and 0„,ax = 0.58 denote the minimum and maximum filling fac- 
tors, respectively, in the compressive strength relation. However, 
the material can exceed both of these values. The power of the 
compressive strength relati on is In(lO) time s the parameter A 
with A - 0.58 (see Gerets hauser et al.ll2010i for a study of this 
parameter), which can be found together with the other mate- 
rial parameters in Table [T] For <p < <p,^in + s, the compressive 
strength relation is continuously extended by the constant func- 
tion 2(0) = 2(0n,in + e) and for (p^^^ < we set 2(0) - oo. The 
tensile strength is given by 

r(0) = -10^'*+' ''**'*Pa , (6) 

and the shear strength is the geometric mean of the compressive 
and tensile strength 

r(0) = ^JnWm ■ (7) 

Hence, the full equation of state for the hydrostatic pressure 
reads 



(2(0) for0+<0 
p(0)= K(cf>'^) for0,T<0<. 
r(0) for0<0- 



(8) 



where 0^ > 0^, and 0^ and 0^ are critical filling factors. The 
value of 0^ marks the transition between elastic and plastic com- 
pression and 0^ defines the transition between elastic and plastic 
tension. The filling factors in-between the critical values rep- 
resent the elastic regime. The quantity 0q is the reference den- 
sity at vanishing external pressure and Kiip) is the filling-factor- 
dependent bulk modulus given by 



(9) 



where y - 4,Ko is the bulk modulus of an uncompressed random 
ballistic deposition (RBD) dust sample with 0rbd = 0.15 (e.g. 
[Blum & Schrapler2004) such that = ^:(0rbd) = 4.5 kPa, and 
the shear modulus is given by /y(0) = 0.5K{(p). 

Plasticity is modelled such that once the filling factor ex- 
ceeds the critical values <p^ or 0^7, the hydrostatic pressure, 
which is computed according to the elastic equation of state 
(second line in Eq. [8]l, is reduced to 2(0) or r(0), respectively. 
For the reduction in pure shear, the deviatoric stress tensor Safi 
is hmited according to the van Mises plasticity given by (e.g. 
iBenz&Asphaugll 19941) 



(10) 



where / = min |^ 7^/372, l], Ji - Sa/sSafi is the second irre- 
ducible invariant of the deviatoric stress tensor, and Y = 7(0) is 
the shear strength. 



Physical Quantity 


Symbol 


Value 


Unit 


Bulk modulus 




4,500 


Pa 


Bulk modulus exponent 


r 


4 




Matrix density 


Ps 


2,000 


kgm 


RBD filling factor 


0RBD 


0.15 




E mean pressure 


Pm 


260 


Pa 


S max. filling factor 


0max 


0.58 




Z min. filling factor 


0min 


0.12 




Z power parameter 


A 


0.58 





Table 1. SiOa dust material parameters used in the porosity 
model. RBD stands for the random ballistic deposition dust sam- 
ple bv lBlum & Schrapled (|2004 and 2 = 2(0) is the compres- 



sive strength. 



2.2. Calibration procedure 

The above porosity model was carefully calibrated and tested in 
IGeretshauser et al! (1201 Ol) wit h the aid of three benchmark ex- 
periments (iGiittler et al. I I2OO9I) : ( 1 ) the compaction of dust by a 
glass sphere, (2) the rebound of a dust aggregate from a solid 
surface, and (3) the fragmentation of a dust aggregate colliding 
with a solid surface. This is the most elaborate material calibra- 
tion which is available for astrophysically relevant macroscopic 
dust aggregates. All benchmark experiments were carried out in 
the laboratory as well as in simulations. 

Because compaction of porous material is a very efficient en- 
ergy dissipation mechanism that allows dust aggregates to stick 
rather than break apart, in our first experiment we tested the com- 
pression properties of the model, ensuring that both the dynamic 
compressive strength and the shear strength relations were cali- 
brated. 

Since rebound is discussed as an obstacle to planetesimal 
formation, the second benchmark experiment tested the bounc- 
ing properties of a dust aggregate. With the aid of this benchmark 
setup, we calibrated the bulk modulus as a characteristic quan- 
tity of the elastic properties. We also confirmed the relation of 
the compressive strength which is a characteristic quantity for 
the energy dissipation during compaction. 

For pre -planetesimal collisions that end in fragmentation, the 
largest fragment and the mass distribution of the fragments is 
highly important for subsequent collisions. For this reason, our 
last calibration experiment tested the correct reproduction of the 
fragmentation properties of the porosity model. The bulk modu- 
lus, which is the quantity that has the strongest influence on the 
fragment mass distribution, was also confirmed. 

After the successful calibration procedure, we are able to 
correctly reproduce the compaction, bouncing, and fragmenta- 
tion behaviour of porous SIOt dust aggregates quantitatively 
with one consistent set of material parameters as shown in Table 
[U These material parameters are used for all simulations dis- 
cussed in the following. 



3. Reproducing sticking, bouncing, and 
fragmentation collision types 

IGiittler et al.l (I2OIOI) compiled a large number of experiments 
with dust aggregates and suggested categorising them into four 
types of sticking, two types of bouncing, and three types of frag- 
mentation. The sticking events were subdivided according to 
sticking mechanisms (hit-and-stick for micron sized monomers, 
surface effects, or deep penetration) and a quantitative crite- 
rion (mass transfer). Growth neutral bouncing with compaction 
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Bl (bouncing with compaction) F2 (erosion) 

Fig. 1. The outcome o f our simulations of pre-planetesimal collisions encompasses all sticking (S), bouncing (B), and fragmentation 
(F) types proposed bv lGiittler et al.l(l2010l) . The initial configuration for each simulation is a sphere with radius rt = 10 cm as resting 
target and a sphere with rp = 6 cm as projectile except for S3 and F2, where - 2 cm. The colour code indicates the filling factor (p. 
Both objects are initially set up with (p = 0.35 except for the F2 case, where 4> = 0.55. The simulations are carried out with different 
impact velocities and the snap shots are taken at different times. The details are given in the text. 



and bouncing with mass transfer were distinguished in terms of 
rebound. Fragmentation was classified according to its degree 
(complete fragmentation or erosion) and whether some sticking 
was involved (fragmentation with mass transfer). 

In this section we show that our code and its underlying cal- 
ibrated porosity model can reproduce these empirical outcome 
types in numerical simulations. Images of the final result of each 
example simulation are presented in Fig.[T]according to this cate- 
gorisation. The end times of the simulations are chosen by visual 
control of the fragment evolution and particle velocities close to 
the thermal equilibrium for each fragment. For each category, an 
exterior view (left) and a cut through the centre of the resulting 
aggregate(s) (right) is shown (except for type Fl). The target is 
always a dust aggregate with radius rt = 10 cm modelled us- 
ing 238,238 SPH particles. The projectiles are modelled using 
1,905 and 51,477 SPH particles for projectile radii rp = 2 cm 
and Tp = 6 cm, respectively, and placed on a cubic lattice with 



edge length 2.6 mm. The collision partners are aligned head on 
with a shift of half a lattice constant into each direction perpen- 
dicular to the collision axis to avoid particle interpenetration. In 
all simulations (except F2), both objects are set up with (p - 0.35 
and the masses of the target, and the 6 cm-, and 2 cm-projectiles 
are 2.93 kg, 0.63 kg, and 0.023 kg, respectively. An overview 
of the initial conditions of the simulations is given in Table |2] 
The colour in Fig. [T] indicates the filling factor (p. We now dis- 
cuss the simulation outcomes in detail. We adopt the notation by 
lGuttleretal.l(l20Tol) . 

SI (hit & stick) is not represented in Fig.[T]as it is not applica- 
ble within the continuous SPH scheme. This category describes 
the outcome of collisions in which microscopic fractal dust ag- 
gregates are involved. In this regime, the continuum limit, which 
is a fundamental assumption of our model, is invalid. We focus 
on collisions of macroscopic dust since for the fractal regime the 
numerical and empirical basis is profound. 
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52 (sticking through surface effects) The example shows the 
result of the impact of an aggregate with radius r^, - 6 cm and 
collision velocity vq = 2.0 ms"'. The image shows the situation 
250 ms after the impact. The filling factor of the outer shell of 
the aggregates remains nearly unchanged as indicated by the ex- 
terior view. In contrast, large parts of the interior are compressed 
to ~ 0.45. The target and the projectile merge to form one ob- 
ject. The initial setup matches that used for the B2 case but has 
a higher impact velocity. As a consequence, the contact area is 
larger and compressed to higher filling factors. This leads to a 
higher tensile strength in this region, which prevents the objects 
from rebounding from each other This situation is referred to as 
sticking through surface effects. A single aggregate remains and 
no fragments are ejected. 

53 (sticking by deep penetration) is found in our simulations 
if the projectile is sufficiently smaller than the target. In the ex- 
ample case, the projectile has a radius of rp = 2 cm and hits the 
target with vo - 10.0 ms"'. The image is taken 388 ms after the 
impact. A small crater is formed on the target and a small num- 
ber of single SPH particles are ejected. This characterises deep 
penetration. The crater is visible in the cross-section. Material is 
compressed to a maximum of (p ~ 0.52 down to a depth of about 
one target radius. 

54 (mass transfer) occurs if the projectile is large and moves 
fast enough to stick but still slow enough to not fully disrupt 
the target. As an example, we choose an impact of a 6 cm ag- 
gregate onto a target with the same filling factor The impact 
velocity is vo — 10.0 ms"' and the image shows the result after 
50 ms. During the collision, the target is highly deformed and 
compressed to filling factors of ^ ~ 0.45 to 0.52. Single SPH 
particles are ejected. Most of the projectile merges, which ac- 
counts for the mass transfer The final mass of the largest frag- 
ment is 3.56kg. 

Bl (bouncing with compaction) is only seen in the case of 
objects of medium and low porosity at low collision speeds. 
Highly porous objects have an extremely low compressive 
strength. As a consequence, low pressures suffice to plastically 
deform the target and the highly porous aggregates do not gain 
enough elastic loading for bouncing. In our example, the projec- 
tile is 6 cm in radius. It impacts at vq = 1.0 ms"'. In Fig.[T] the 
final state at 500 ms is shown. In the collision, both objects are 
flattened at the impact site and the elastic loading is sufficient to 
make them rebound after the impact. Part of the interior of the 
target and projectile is compressed to ~ 0.45. We note that this 
is exactly the same setup as in the S2 case but with a lower im- 
pact velocity. The bouncing event is a result of a smaller contact 
area and a lower tensile strength in this region, which is indicated 
by the lower filling factor As a consequence, the two objects re- 
bound. A close investigation reveals that in nearly every simu- 
lation with bouncing, a small amount of material is transferred 
between the projectile and target. Thus, it is very unlikely that 
pure growth neutral bouncing without any mass transfer exists. 
In this particular simulation, material of ~ 0.54 g is transferred 
from the projectile to the target. Hence, the result could also be 
classified into the next category. 

B2 (bouncing with mass transfer) occurs in simulations 
where the compressive strength is sufficiently large to allow for 
elastic loading and consequential bouncing, and where the ten- 
sile strength is small enough for the impactor to rip out a small 
fraction of mass from the target. This happens at very low im- 
pact speeds. The final state at 1.7 s after the collision is the result 
of a collision with vq = 0.2 ms"'. The projectile has r^ = 6cm. 
During the impact, a small region of target and projectile is com- 



pressed to ~ 0.45. In the consequential rebound, this region 
sticks to the projectile and is ripped out of the target. The re- 
maining crater on the target can be seen in the cross-section. In 
this collision, 29.5 g are transferred from the target to the projec- 
tile. 

Fl (fragmentation) is generally the outcome of collisions 
with high impact velocities. Highly porous objects effectively 
dissipate energy by deformation because of their low compres- 
sive strength but they are also easy to disrupt because of their low 
tensile strength. In contrast, objects of low porosity have high 
tensile strengths, but lack the ability to dissipate large amounts of 
kinetic energy. The degree of fragmentation is therefore strongly 
porosity dependent. Our example shows the result of a collision 
between objects, where the projectile features r^ - 6 cm. The 
impact speed is vq = 17.5ms"'. During the collision, both the 
target and the projectile completely shatter and the result of this 
collision (shown only in exterior view 800 ms after the impact) 
is a continuous fragment distribution, whose masses range from 
285 g down to single SPH particles of 12 mg. The fragments con- 
sist of target and projectile material that are indistinguishably 
combined together The filling factors of the fragment distribu- 
tion are ~ 0.54, which is close to the maximum filling factor 
of 0.58 (seeTableffli. 

F2 (erosion) is observed particularly for high filling factors 
and small projectile radii. The initial setup for the example case 
involves two objects with = 0.55. The target and projectile 
masses were 4.60 kg and 0.037 kg, respectively. The projectile 
radius is r^ - 2 cm and the impact velocity vq = 20.0 ms"'. 
During the intrusion of the projectile, fragments consisting of a 
small number of SPH particles are ejected opposite the impact 
site. The figure shows the situation 438 ms after the collision, 
where small fragments and SPH particles are ejected from the 
crater In this collision, 63.2 g of material are eroded from the 
target. The projectile intrudes to the centre of the target and the 
intrusion channel is clearly visible in the cross-section. Because 
of the high compressive strength elastic deformation also takes 
place on the target. The figure shows a snap shot of a decom- 
pression wave of lower filling factor, which appears as a vertical 
ring on the target. 

The category F3 (fragmentation with mass transfer) con- 
tains conceptual difficulties. The first reason is that in the sim- 
ulation outcomes of the Fl and F2 categories, some mass is 
always transferred to the target. Secondly, the demarcation be- 
tween sticking with mass transfer (S4) and F3 is unclear and a 
continuous transition between these categories is expected. This 
problem can be solved by our new model proposed below. Our 
simulations incorporate category F3 by reproducing the S4 and 
Fl types. 

In conclusion, our code is not only capable of quantitatively 
reproducing sticking, bou ncing, and fragmentation in general 
dGeretshauser et al]|20I0h but also of correctly simulating the 
sub-types for each macroscopic collision outcome. However, we 
also experienced difficulties in cla ssifying the results o f our sim- 
ulations according to the model of lGiittler et al.l (|20I0|) . This mo- 
tivates our wish to find a new approach. 



4. A new model for mapping collision outcomes 

We introduce a new model to classify the outcome of pre- 
planetesimal collisions. We first describe its structure. We then 
show that all sticking, bouncing, and fragmentation events of the 
previous section can be incorporated into the new model. 
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Typevo [ms '] 


r, [cm] 


nit [kg] 


Tp [cm] 


'Wp [kg] 


<^t/p 


S2 2.0 


10 


2.93 


6 


0.63 


0.35 


S3 10.0 


10 


2.93 


2 


0.023 


0.35 


S4 10.0 


10 


2.93 


6 


0.63 


0.35 


Bl 1.0 


10 


2.93 


6 


0.63 


0.35 


B2 0.2 


10 


2.93 


6 


0.63 


0.35 


Fl 17.5 


10 


2.93 


6 


0.63 


0.35 


F2 20.0 


10 


4.60 


2 


0.037 


0.55 



Table 2. Initial parameters of the simulations shown in Fig. [T] 
The quantities are: collision velocity vq, target radius rt and mass 
nit, projectile radius rp and mass nip, and the filling factor of the 
collision partners <pt/p. 



4.1. The motivation beliind a new classification sclieme 

The model o f lGiittleret all ( 120101) was developed based on par- 
ticular laboratory setups. In several of these, one of the collision 
partners was not a dust aggregate but a solid surface or a glass 
bead, which itself cannot fragment. Therefore, applying this cat- 
egorisation to pre-planetesimal dust collisions, which here are 
solely carried out as simulations between dust aggregates, leads 
to some difficulties. 

In our simulations, the collisions of S3 (sticking by pen- 
etration) always produce ejected dust. Therefore, it is unclear 
whether they should be classified as S3 or rather S4 (mass trans- 
fer). On the other hand, S2 (sticking through surface effects) in- 
volves also a mass transfer but without the production of frag- 
ments. It can also be expected that there is a continuous transi- 
tion from S2 to S3 with increasing impact velocity. The demar- 
cation between the sticking types is also conceptually difficult to 
make. The types SI to S3 are distinguished on the basis of the 
sticking mechanism, whereas for S4 the criterion is an increase 
in target mass and the generation of some fragments. In the con- 
text of pre-planetesimal growth, the exact sticking mechanism is 
of minor importance. The distinction between growing and dis- 
ruptive events is simply given by comparing the largest fragment 
before and after the collision. 

Applying the bouncing categories to collision data also 
causes some difficulties. Pure bouncing with compaction (Bl) 
events are never seen in our simulations. Analysing the final 
mass of both collision partners after the impact reveals that some 
mass is always transferred either from the projectile to the tar- 
get or vice versa. In the first case, growth occurs and the re- 
sult should be categorised under S4 with a growing target, but 
with one remaining fragment instead of a fragment distribution. 
If mass is transferred to the projectile, then the event belongs 
to B2 (bouncing with mass transfer). However, since the largest 
object is losing mass, B2 is also a type of erosion (F2), where 
instead of a fragment distribution only one large fragment (the 
enlarged projectile) is present. Since bouncing m ay be an ob- 
stacle to planetesimal growth (IZsom et al ] 1201 01) . the effect of 
bouncing has to be included in any colUsion map. However, in- 
stead of distinguishing between two types of bouncing, we find 
that it is sufficient to distinguish between the mass of the largest 
and the second largest fragment. 

Analysing the numerical results for disruptive events, it is 
also hard to distinguish between types Fl (fragmentation) and 
F3 (fragmentation with mass transfer). As already mentioned 
above, in any fragmenting collision, mass is transferred from the 
projectile to the target and the fragments consist of both pro- 
jectile and target material. In eroding events (F2), mass is also 
transferred to the target and a fragment distribution is produced. 



Once again the continuous transition between the three fragmen- 
tation types and between sticking and fragmentation events can- 
not adequately be described by the given model. In our view, it 
is sufficient to characterise the outcome of disruptive events in 
terms of the size of the largest fragment (which for F2 is much 
larger than the other fragments) and a power-law distribution of 
the remaining fragments. Thus, for correct mapping of any com- 
bination of sticking, bouncing, and fragmentation, the following 
objects have to be considered: both a largest and a second largest 
fragment as well as a power-law fragment distribution. 

This discussion shows that the classification of different 
sticking, bouncing, and fragmentation events is not sufficient. 
This is because many intermediate events exist. For example, 
a bouncing event can involve some sticking and fragmentation, 
whereas a stic king event can also involve some fragmentation. 
The model by iGiittler et al.l ( 1201 Ol) tries to capture these inter- 
mediate events. However, it is based on physical mechanisms 
rather than collision outcome. Owing to its discrete approach, 
the model generates some unnecessary complexity and cannot 
easily describe transitions between categories. To improve on 
this, we propose a model that is based solely on quantitative as- 
pects. It distinguishes between four types of fragment popula- 
tions characterised by continuous quantities such as their mass, 
filling factors, velocities, and size. In this way, we are able to 
map all types of sticking, bouncing, and fragmentation events 
but also model continuous transitions between these types. 

4.2. Fragment populations 

The fundamental idea of the new classification is that any out- 
come of a collision can adequately be modelled by distinguish- 
ing four kinds of fragment "populations". This is illustrated in 
Fig. |2] Example input parameters for the collision include the 
masses of the target (label "t") nit and projectile (label "p") Wp, 
their filling factors 0t and 0p, their rotational energies Eyox.t and 
/ii-ot.p, their impact parameter b and, most importantly, their col- 
lision velocity vq (see Fig.|2] left). Depending on these parame- 
ters, the collision between two bodies produces a well-defined 
outcome of populations. To be able to map sticking, bounc- 
ing, and fragmentation events, it is suitable to sort this distri- 
bution according to the mass of the fragments in the respective 
population. For these, output parameters such as mass, filling 
factor, rotational energy, and velocity are determined as exact 
or averaged values, or as distributions depending on the pop- 
ulation. Specifically, the new scheme consists of the following 
four types: (1) largest fragment, (2) second largest fragment, (3) 
power-law population, and (4) sub-resolution population. 

(a) Largest fragment: Comparing the mass of the largest 
fragment before and after a collision distinguishes between posi- 
tive and negative growth. The characteristic output quantities are 
supplied. 

(b) Second largest fragment: This enables bouncing events to 
be mapped. In pure bouncing events, the largest and the second 
largest fragment will be the only members of the fragment popu- 
lation. The quantities describing this fragment are supplied. The 
second largest fragment only exists if it consists of more than a 
single SPH particle. 

(c) Power-law population: in disruptive collision events in 
particular, the mass distribution of the fragments can be mod- 
elled by a power-law to describe the fragment distribution. 
Instead of describing each fragment of the power-law popula- 
tion with quantities as done previously for the largest and sec- 
ond largest fragment, the number of parameters is reduced sig- 
nificantly by utilising distribution functions. 
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Fig. 2. Illustration of the four-population model. The left side displays the situation before the collision: two pre-planetesimals 
collide with impact velocity vq and impact parameter b. The target and projectile are characterised by their defining quantities, 
e.g., mass m, filling factor (p, and rotational energy Ei-ot- Depending on these parameters, an outcome population is generated in 
the collision (middle). The four-population model distinguishes the fragments according to their mass. The population classes are: 
(a) the largest fragment, (b) the second largest fragment, (c) the power-law population, and (d) the sub-resolution population. This 
categorisation is sufficient to describe the outcome types of Fig.[T]as well as intermediate outcomes. 



(d) Sub-resolution population: This population is introduced 
for numerical reasons. Fragments of the sub-resolution popu- 
lation consist only of a single SPH particle. The existence of 
this population represents the resolution limit of the simulations 
and therefore represents only an upper limit to the smallest dust 
fragments, which are produced in a collision. This population is 
also convenient to control sufficient resolution. The SPH numer- 
ical scheme is capable of simulating objects of metre size and 
more. However, because of limited computational resources, sin- 
gle SPH particles then represent objects of centimetre size and 
larger. The sub-resolution population keeps track of these ob- 
jects. It is convenient to describe this population with averaged 
characteristic quantities. We note that also in laboratory experi- 
ments a sub-resolution population exists, for which the fragment 
distribution is determined by high speed cameras. 

With this approach, we considerably reduce the complexity 
of the sticking, bouncing, and fragmentation classification from 
nine types of collision outcomes to four fragment populations. 
Moreover, our model is based solely on continuous quantities 
characterising the collision outcome and not on physical mech- 
anisms. This allows us to model any mixed types of collisions 
and transitions between growth and disruption with the neces- 
sary accuracy. 

Compared to fragmentation data t hat are mapped to a power- 
law distribution alone (e.g. Mathis et al.l 19771 iDavis & RvanI 
ll99nllBlum & Munchlll993irGuttler et al.ll2010l) . the new model 
is also more accurate. In some of our simulations, a fragment dis- 
tribution contains one large fragment and many small ones. The 
small ones can be accurately modelled with a power-law mass 
distribution. However, the mass of the largest fragment often 
does not match the power law. This is particularly true for graz- 
ing collisions between highly porous aggregates, where the fill- 
ing factor of the largest remnant differs significantly from those 
of the other fragments. Since this has a major impact on sub- 
sequent collisions, we treat the largest fragment and its filling 
factor separately. Furthermore, power-law distributions cannot 
map bouncing collisions. For these reasons, we distinguish the 



largest and second largest fragment from the rest of the fragment 
distribution. 

With the four given fragment populations and their defining 
properties, we present a closed model that is capable of mod- 
elling any collision outcome of pre-planetesimals with minimum 
complexity but with the accuracy necessary to model the dust ag- 
gregation in global coagula tion models. In t he next section, we 
show how the model bv iGiittler et al.l (1201 Oh can be represented 
in our new model. 



4.3. Mapping sticking, bouncing, and fragmentation to tine 
new modei 

The key idea behind mapping the sticking, bouncing, and frag- 
mentation sub-types to the four-population model is to describe 
the collision outcome by the characteristic quantities of each of 
the four populations. If one of the populations does not exist, its 
characteristic quantity is simply set to zero. In principle, a num- 
ber of quantities can be used to determine the population. To 
illustrate the use of this model, we select the characteristic quan- 
tities mass m and filling factor <p as e xamples. Table [3] d emon- 
strates the successful mapping of the Guttler et al.l (1201 Ol) types 
into our quantitative four-population model. The indices "t", "p", 
"1", "2", "pw", and "sr" denote target, projectile, largest, second 
largest, power-law, and sub-resolution quantities, respectively. 
For the illustration, we assume that the mass of the target nii is 
greater than or equal to the projectile mass nip, i.e. nit ^ nip. We 
find a filling factor increase in all collisions if the initial filling 
factor is smaller than the maximum compaction at i^^ax - 0.58 
(see Fig.Qfor the largest fragment). 

For pure sticking (SI to S3 types) only one fragment in the 
final population exists, which is identified as the largest frag- 
ment. Consequently, this fragment contains the mass of the to- 
tal system, i.e. nii = + nip. Within the four-population ap- 
proach, the SI to S3 sub-types can be combined in terms of the 
mass and filling factor In general, the filling factor of the tar- 
get <pt is increased during the impact such that the filling factor 
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Table 3. Mapping sticking, bouncing, and fragmentation. In this table, we illustrate how the collision outcome types of Sec.|3]can 
be mapped to the four-population model using the characteristic quantities mass m and filling factor (p as examples. The masses 
of the largest and second largest fragment, the power-law population, and the sub-resolution population are given by m\, 1112, nip^,, 
and msi, respectively. We assume for the target mass nit and projectile mass nip before the collision nit ^ '«p- The filling factors 
are represented by 0t for the target, (pp for the projectile, (pi for the largest fragment, and 02 for the second largest fragment. We 
also give examples of the transition criteria between the types of Sec. [3] A list of symbols can be found in the appendix and further 
explanations are given in the text. 



outcome type 


ITl] 


4>i 


;;72 


02 








hit & stick (SI) 


nit + "'p 


4>i 
















sticking through surface effects (S2) 


rtit + OTp 


01 >0t 
















sticking by penetration (S3) 


m, + OTp 


<Pi ><Pt 
















mass transfer (S4) 


mi > rut 


01 >0t 


m2 < OTp 


02 > 


Pv 


Wlpw > 


ms, > 


bouncing with compaction (Bl) 


mi = mt 


01 >0t 


m2 = nip 


02 > 


Pv 








bouncing with mass transfer (B2) 


mi mt 


01 >01 


m2 t nip 


02 > 


P9 








fragmentation (Fl) 


mi < mt 


01 >01 


m2 > 


02 > 


Pp 


'Wpw > 


m^r > 


erosion (F2) 


mi < mt 


01 >0t 


m2 > 


02 > 


Pp 


Wlpw > 


ms,- > 


fragmentation with mass transfer (F3) 


mi > mt 


01 >0t 


m2 > 


02 > 


Pp 


JWpw > 


ms,. > 


sticking - fragmentation transition 


> mt — > mi < mt 


01 


m2 > 


02 




mp„ > 


msr > 


bouncing - fragmentation transition 


mi < mt 


01 


m2 > 


02 




tnpw = — » OTpw > 


OTsi = — > OTsr > 



of the largest fragment <pi > 4>i- For sticking with mass trans- 
fer (S4), a range of fragments exists. Since sticking is identified 
with growth of the largest fragment, it is mi > mt. If there are 
only fragments below the resolution limit, then the mass of the 
second largest fragment mi and the mass of the power-law popu- 
lation mpw vanish. Otherwise there are a range of fragments with 
non-zero masses. 

Bouncing is characterised by two fragments in the final 
fragment distribution. For bouncing with compaction (Bl), the 
masses are unaltered, i.e., ni\ - nit and ni2 - nip, but the filling 
factors are increased such that 0i > 0t and 02 > 0p- For the case 
of bouncing with mass transfer (B2), the filling factors are also 
increased but the masses are altered such that m\ + ni2 - nit + nip 
and m\ + nii well as m2 nip. By definition, no fragments 
are produced in the bouncing events such that mpw - msr - 0. 
For B2, the filling factors are also increased, i.e., 0i > 0t and 

02 > 0p- 

The fragmentation events (Fl to F2 types) are characterised 
by the mass of the largest fragment being lower than the highest 
mass before the collision, i.e. nii < mt. The distinction between 
Fl and F2 in terms of the collision outcome is unclear. In both 
cases, a power-law as well as a sub-resolution population is gen- 
erated, i.e. mpw > and msr > 0. It seems that for F2, the mass of 
the target is only re duced by a small amount such that mi < mt 
(lGuttleretal.l2010h . In this reference, F3 is also characterised by 
the target (solid plate) gaining mass, i.e. mi > mt. This is incon- 
sistent with the idea that fragmentation is equivalent to nii < nit. 
The mass of the second largest fragment has some value ni2 > 
and total masses of the power-law and sub-resolution popula- 
tions are mp„ > and ni^r > 0. 

The transition between sticking and fragmentation can con- 
sequently be characterised by the transition mi > mt — > mi < mt, 
since for the S4 type power-law and sub-resolution population 
are already present. In contrast, the transition from bouncing 
(which also includes the change of target and projectile masses 
by mass transfer in B2) can be defined by the appearance of a 
power-law and sub-resolution population, i.e. mp^ = — > mp„ > 
and msr = — > msr > 0, together with a non-growing target 
mi < mt. 

In this section, we have shown that the four-population 
model for collision outcomes is capable of encompassing all the 
sticking, bouncing, and fragmentation sub-types. Furthermore, 



Collision 


Broad 


Velocity [ms"'] 


categorisation 


0.1 


bouncing 


0.2 


bouncing 


1.0 


bouncing 


1.25 


sticking 


1.5 


sticking 


1.75 


sticking 


2.0 


sticking 


2.2 


sticking 


2.5 


sticking 


4.0 


sticking 


5.0 


sticking 


6.0 


sticking 


7.5 


sticking 


10.5 


sticking 


11.5 


fragmentation 


12.5 


fragmentation 


15.0 


fragmentation 


17.5 


fragmentation 


25.0 


fragmentation 


27.5 


fragmentation 



Table 4. Collision velocities used in the simulations. The tar- 
get and projectile are homogeneous and feature an intermediate 
porosity (0 = 0.35). These simulations have been carried out 
to show how the proposed model can quantitatively demonstrate 
the results of pre-planetesimal collisions. Each simulation is as- 
signed a broad categorisation so that the transition regions be- 
tween the categories can be investigated. 



we have demonstrated that the transition between these types 
can be modelled continuously by utilising the masses of each of 
the four populations. 

5. Applying the new model to simulation data 

We apply the four-population model to a study of collisions 
of medium porosity pre-planetesimals with different velocities 
to demonstrate that the outcomes of pre-planetesimal collisions 
can be described quantitatively using the model. We caiTy out 
24 simulations of head-on collisions involving spherical dust 
aggregates with target and projectile radii of rt = 10 cm and 
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-0.5 0.5 

log collision velocity [ms" 



1.5 



Fig. 3. Cumulative plot of the contributions to the total mass in 
the system by the largest (grey slashed), second largest (solid 
grey), power-law (black slashed), and sub-resolution (solid 
black) populations after the aggregate collisions at various ve- 
locities. At low velocities, where bouncing occurs, the mass is 
wholly within the first and second largest fragments, while at 
intermediate velocities, where sticking occurs, the contribution 
to the mass is in the largest fragments. At high velocities, the 
mass contribution from the power-law population becomes sig- 
nificant. 



Tp - 6 cm modelled using 238,238 and 51,477 SPH particles, 
respectively. We choose the initial filling factor of the aggre- 
gates to be - 0.35 resulting in target and projectile masses 
of 2.93 and 0.63 kg, respectively. We vary the collision velocity 
between 0.1 and 27.5 ms"'. Table |4] summarises the simulations 
carried out. We broadly categorise the simulations by eye into 
bouncing, sticking, ov fragmentation so that a comparison can 
be made between the results of our model and the categorisa- 
tions used previously. We note that these broad categorisations 
are only employed to show where the boundaries between these 
regimes lie and are not used for the quantitative approach of our 
model. 

As discussed in Sec. 14.21 the outcome of pre-planetesimal 
collisions can be quantitatively described by a number of pa- 
rameters such as the final mass, size, energy, velocity, porosity, 
and rotation and each of these parameters can be used to describe 
the four populations. In this preliminary, study we focus on the 
final mass, porosity, and energy of the different populations in 
our model. 

Figure [3] shows how the mass distribution amongst the four 
populations varies with collision velocity. At low collision ve- 
locities, the contact energy and thus the tensile strength is too 
low to hold the objects together This results in bouncing caus- 
ing the overall mass distribution being similar to the initial dis- 
tribution. The collision with 0.2 ms"' shows evidence of mass 
transfer resulting in the largest fragment being smaller than the 
target. As the collision velocity increases, the energy is dissi- 
pated by plastic deformation resulting in sticking such that all 



the mass is stored in a single object. At even higher collision ve- 
locities in which fragmentation is seen, the mass stored in the 
largest fragment is reduced but is increasingly present in the 
power-law population. Figure |3] also shows how the mass of the 
largest fragment varies with collision velocity. It can be seen that 
three distinct regions exist, which justify the broad classifica- 
tion of Table H] However, there do not appear to be any distinc- 
tions within these regions which justify an additional division 
into sticking, bouncing, and fragmentation sub-types. A sharp 
transition between the bouncing and sticking regions exists at 
~ 1 ms"'. Above 10.5 ms"', pure sticking no longer occurs and 
the mass of the largest fragment begins to decrease as the transi- 
tion between sticking and fragmentation takes place. As a result 
of very high collision velocities, the mass contribution of the 
sub-resolution population dominates over the largest and second 
largest fragment which is an indicative of violent disruption. We 
note that the simulation with collision velocity 11.5 ms"' still re- 
sults in collision growth but above this velocity the mass gain of 
the largest fragment is negative. 

Figure [3] clearly shows that at velocities > 11.5 ms"', pure 
sticking no longer occurs and a transition into the fragmenta- 
tion regime begins. For the simulation with a collision velocity 
of 11.5 ms"' (see Fig.|4] top, for an illustration of the resulting 
fragments), the mass of the l argest fragment clea rly increases 
but according to the model bv lGiittler et al.l (l20I0l) . it is unclear 
whether this simulation would be classed as sticking with mass 
transfer (S4) or fragmentation with mass transfer (F3). Using 
our model, such an intermediate region can be described quanti- 
tatively. At ~ 12.5 ms"' (see Fig.H] bottom, for an illustration), 
the mass contribution from the power-law population becomes 
as significant as the contribution from the largest fragment. In 
this region, the target fragments into a small number of large 
pieces. We expect that the transition from the sticking to frag- 
mentation regime, described in Section l43] occurs between 11.5 
and 12.5 ms"'. 

Using the model, it is clear that the mass of the largest frag- 
ment decreases sharply for collision velocities > 1 1 .5 ms"', thus 
this is the region where fragmentation takes place. We expect 
that this fragmentation boundary will change depending on a 
numb er of other parameters such as impact parameter and mass 
ratio dBlum & Miinchiri993h . porosity, and the rotation of the 
dust aggregates. Further investigation into the factors that deter- 
mine the fragmentation boundary is crucial to ultimately under- 
stand under what conditions fragments may grow to planetesi- 
mal sizes. It is important to note that to decrease the computa- 
tional expense, the SPH particles that move out of a radius of 
20 m are removed from the simulation. Therefore, the total mass 
at the end of the higher velocity simulations is smaller than the 
initial mass. However, the total mass removed from the simula- 
tion is ~ 2% in the highest velocity simulation. 

Figure|5]shows how the mass of the second largest fragment 
varies. For low collision velocities, it is evident that bouncing 
occurs and the projectile, whose mass is indicated by the dashed 
line in Fig. |5j appears as the second largest fragment with un- 
altered or slightly increased mass in the final distribution. For 
higher collision velocities, the projectile sticks to the target and 
no second largest fragment results, which is indicated by van- 
ishing masses in this regime. As the velocity further increases, 
the resulting mass of the second largest fragment increases again 
but is smaller than the projectile mass, which indicates fragmen- 
tation. The increase is followed by a sharp decrease in the mass 
of the largest fragment. These simulations, which are reasonably 
close to the sticking-fragmentation boundary not only cause the 
two aggregates to fragment, but also cause mass to be transferred 
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Fig. 4. Outcome of the dust collision simulation with 11.5 ms ' 
(top) and 12.5 ms"' (bottom) at a time, t=0.8 s. The collision 
with 1 1 .5 ms"' results in the growth of the target but also a large 
chunk does break off such that pure sticking does not occur The 
situation is adequately mapped with the appearance of a second 
largest fragment. For 12.5 ms"' the target breaks apart and the 
power-law population becomes as significant as the largest frag- 
ment. The transition from sticking to fragmentation occurs in 
between the presented collision velocities. 



between the target and projectile. At higher velocities still, the 
impact is sufficiently violent for the second largest fragment to 
also decrease sufficiently. 

Figure [3] shows that as the collision velocity increases, the 
contribution to the mass from the power-law population also in- 
creases. As an outcome of pre-planetesimal collisions, the cu- 
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Fig. 5. Mass of the second largest fragments. The broad categori- 
sation into bouncing (open triangles), sticking (closed triangles), 
and fragmentation (open squares) can be identified by masses 
of the second largest fragments being approximately the projec- 
tile mass (dashed line), nearly vanishing, and smaller than the 
projectile mass, respectively. For the fragmenting collisions, the 
mass of the second largest fragment initially increases with col- 
lision velocity. However, at higher collision velocities, the frag- 
ment mass decreases because of the collision nature of the sim- 
ulations. 

mulative mass distribution of the fragments is often described 
by a power-law 



mcumim) 



Zr"" I nif Y 
nif — I n{m)mdm-\ , (11) 



where n(m) dm is the number of fragments within the mass range 
[m, m+dm], nif is the fragment mass, and yUp„ is the most massive 
member of the power-law population. The quantities m, nif, and 
fipv, are normalised by the total mass of the power-law population 
and k is the power-law index. 

Figure |6] shows the cumulative mass distribution mcum('Mf) 
against the fragment mass normalised by nip^. For higher frag- 
ment masses nif, a power-law fit can be obtained that describes 
the mass distribution of the power-law population. As the ve- 
locity increases, the slope of the power-law decreases. This is 
because at higher collision velocities, the destruction of the dust 
aggregates is more violent and a larger number of smaller frag- 
ments results. At low fragment masses, it is currently unclear 
whether the deviation from the power-law distribution is phys- 
ical, or whether it is a result of a low number of SPH particles 
(< 100) per fragment. We note that this deviati on is also seen in 
the experimental results of iGiittler et al.l (l2010l) . 

We now consider how much compaction takes place in each 
of the simulations by considering the resulting filling factors of 
the various populations. Figure |7]shows the average filling factor 
of the largest fragment, (pi, compared to the initial filling factor 
As the velocity increases, the compaction also increases, causing 
the filling factor to increase with velocity. This forms a smooth 
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Fig. 6. Cumulative mass distribution (given by equation fTTT i 
against fragment mass of the power-law populations resulting 
from the simulations with 17.5 (open triangles), 25.0 (closed tri- 
angles), and 27.5 ms"' (open squares). Both axes are normalised 
by TOpw. For high fragment masses, the gradient of the slope de- 
creases with increasing velocity. For low masses, it remains un- 
clear whether the deviation from the power-law is a physical or 
numerical artefact. 
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Fig. 7. Final average filling factor of the largest fragment in 
each of the simulations. The simulations have been broadly cat- 
egorised into bouncing (open triangles), sticking (solid trian- 
gles) and fragmentation (open squares). The filling factor varies 
smoothly as the collision velocity is increased. 
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Fig. 8. Final average filling factor (black line) of the power-law 
populations in each of the simulations. We also plot the max- 
imum (upper grey line) and minimum (lower grey line) filling 
factors for each simulation, as well as the initial filling factor 
(dashed line). As the collision velocity increases, the range of 
filling factors for any one simulation increases. 



curve between the initial filling factor of 0.35 and the maximum 
filling factor of ~ 0.58 for SiOo- The curve resembles the Fermi- 
Dirac shape of the compressive strength relation (Eq. |5), which 
directly links the dynamic pressure to the obtained filling fac- 
tor. We note that the appearance of fragmenting events coincides 
with the largest fragment reaching a filling factor close to the 
maximum filling factor. This suggests that fragmentation sets in 
when the material is maximally compressed. Furthermore, we 
stress that the transition of the filling factor is smooth right from 
the low velocity (bouncing) collisions, through to the medium 
velocity (sticking) collisions to the high velocity (fragmentation) 
collisions. For simulations with higher initial filling factors, we 
expect there to be a direct transition between the bouncing and 
fragmentation regimes, as discussed in Section |43| 

The transition is particularly important because for planet 
formation, we are primarily concerned with the evolution of the 
largest fragment. This demonstrates that our model can quanti- 
tatively capture the results of the collisions over a wide region 
of velocity parameter space. Our preliminary results for the sec- 
ond largest fragment suggest that only a small amount of com- 
paction takes place for the second largest fragment in low veloc- 
ity bouncing simulations such that the final filling factor is very 
close to the initial filling factor, while for high velocity fragmen- 
tation simulations, a large amount of compaction takes place re- 
sulting in filling factors close to the maximum value of 0.58. 
Additional investigations will be carried out in a future study. 

Figure [8] shows the average, maximum, and minimum fill- 
ing factors of the power-law population. At higher velocities 
when more violent fragmentation occurs, the spread in the fill- 
ing factor is greater than at lower velocities. This suggests that 
at lower velocities, fragments may be chipped off without too 
much change to their compaction, whereas in the more violent 
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Fig. 9. Filling factor distribution for the simulations with 17.5 
(solid line), 25.0 (short dashed line) and 27.5 ms"' (long dashed 
line). The size of a filling factor bin is 0.01. As the collision ve- 
locity increases, the number of fragments increases. For higher 
collision velocities, the majority of fragments have lower filling 
factors than the original value (dotted line). 



cases with higher collision velocities, the fragments may become 
compacted before being chipped off or rupture because the plas- 
tic flow occurs where parts of the dust are stretched before be- 
ing ripped off. Figure |9] shows the filling factor distribution of 
the power-law populations resulting from the simulations with 
17.5, 25.0, and 27.5 ms"'. As the collision velocity increases, the 
number of fragments increases. In particular, most of the parti- 
cles have filling factors smaller than the initial value of 0.35. 
This suggests that during the collision, fragments are ripped off 
instead of being compressed before breaking apart. 

Finally, we explore the contributions to the total energy (nor- 
malised by the initial energies) from the different populations. 
The total energy is the sum of the energy stored in transla- 
tion, rotation, and vibration. Figure [10] shows that at low ve- 
locities the major energy contribution comes from the largest 
fragment, while at higher velocities the contribution from the 
power-law population becomes somewhat equally significant (at 
~ 12.5 ms"') and even more so at even higher velocities. At 
low velocities, the energy is stored in elastic deformation and re- 
leased into kinetic energy again so that a smaller fraction of the 
initial energy is dissipated. As the velocity increases plastic de- 
formation occurs causing energy dissipation. At even higher ve- 
locities, though energy is still dissipated (since bonds are broken 
in the fragmentation process), the dissipated energy as a frac- 
tion of the initial energy due to bonds being broken is not as 
large as the energy dissipation due to plastic deformation. For 
very high impact velocities, a considerable amount of energy is 
stored in the sub-resolution population. Figure [TT]shows the ab- 
solute energy dissipated at each collision velocity as well as the 
initial energy. The total amount of energy dissipated increases 
smoothly with increasing velocity, irrespective of the bouncing, 
sticking, or fragmentation regime. We note that in these exam- 



Fig. 10. Cumulative plot of the contributions to the total energy 
in the system by the largest (grey slashed), second largest (solid 
grey), power-law (black slashed), and sub-resolution (solid 
black) populations after the aggregate collisions at various veloc- 
ities. At low velocities, the energy is primarily from the largest 
fragment, whereas at higher collision velocities, where more 
fragments form, the contribution from the power-law population 
becomes more significant. 
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Fig. 11. Initial (dotted line) and dissipated (solid line) energies 
against the collision velocity. The amount of energy dissipated 
changes smoothly with collision velocity, regardless of whether 
the simulations involves bouncing (open triangles), sticking 
(closed triangles), or fragmentation (open triangles). 
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pies, the final rotational energies are not very significant because 
the collisions are head-on. The difference between the rotational 
and translational energies will become more important in simu- 
lations with non-zero impact parameters. 

We have presented here the results of a small number of 
simulations, which show that transitions between the bouncing, 
sticking, and fragmentation regions occur and that the results of 
pre-planetesimal collisions can be described quantitatively. Our 
study could be extended to explore a vast parameter space of in- 
puts. Our model has been shown not only to encompass previous 
qualitative models (SecUJi but also provide an alternative model 
that focuses on the quantitative outcome of pre-planetesimal col- 
lisions. Furthermore, the model can provide an accurate input 
into global models, which may then carry out detailed calcula- 
tions of pre-planetesimal collisions. In addition, this model pro- 
vides a powerful aid to determine the region of parameter space 
that allows growth of pre-planetesimals to occur 

6. Discussion and outloolc 

We have demonstrated that our SPH code, extended and cali- 
brated for the simulation of porous SiOa dust, is capable of re- 
producing all sticking, bouncing, and fragmentation types that 
appear in collision experiments with macroscopic porous dust 
aggregates. In addition to the quantitativ ely correct simula- 
tion of laboratory bench mark experiments (iGiittler et al.ll2009l 
iGeretshauser et al.l20I0l) . this consolidates the validity of the ap- 
plied porosity model and shows its readiness for the application 
in the field of investigating pre-planetesimal collisions. Since the 
continuum approach of the SPH method does not place an upper 
bound on the aggregate size, collision data of a parameter space 
that have been inaccessible so far in laboratory experiments can 
now be provided for further use in global coagulation models. 

For this transfer, a suitable mapping of collision data has 
to be chosen, which is accurate enough to quantitatively cap- 
ture the most important features of any combination of sticking, 
bouncing, and fragmentation and simple enough to be imple- 
mentable in global coagulation models. We have attempted to 
map o ur simulation data to the categorisation by Guttler et id] 
(120 lOl) . which was the most elaborate collision model previously 
available. On the one hand, the distinction between four stick- 
ing, two bouncing, and three fragmentation types introduces a 
considerable amount of complexity caused by distinguishing be- 
tween a mixture of qualitative and quantitative attributes and by 
adhering to a distinct classification of these events. On the other 
hand, we have found a collision outcome that could not clearly 
be attributed to one of the proposed categories. 

Because of this ambiguity, we have proposed a new model, 
which is based on quantitative aspects. For this purpose, we have 
divided the set of fragments of a collision into four populations: 
the largest and second largest fragment are described by distinct 
values for the characteristic quantities of mass, filling factor, and 
kinetic energy to name only a few. The power-law population 
is described by distributions and the sub-resolution population 
by average values of the characteristic quantities. The largest 
fragment indicates growth or erosion, the second largest frag- 
ment accounts for bouncing, the power-law population quan- 
titatively describes the amount of fragmentation, and the sub- 
resolution population gives an upper limit for smaller fragments, 
which are not captured because of insufficient resolution. Since 
the SPH code is not restricted to small aggregate sizes, the im- 
portance of the sub-resolution population becomes significant 
for aggregate collisions between objects of approximately me- 
tre size and more. For growth models that rely on the sweeping 



up of small particles (iTeiser & Wurmll2009l) . the sub-resolution 
population also plays an important role. We have demonstrated 
that this model in g eneral encompasses the model proposed by 
IGiittler et al.l (l2OI0l) but is also capable of capturing intermediate 
events. 

Finally we have applied the new model to map data for 
head-on collisions of aggregates with intermediate {cp = 0.35) 
porosity and varying impact velocity. We have shown that the 
broad sticking, bouncing, and fragmentation categorisation can 
still be found in the four-population model. But in addition, we 
have also shown that continuous transitions in the variation in 
the mass of the largest fragment, the filling factor, and the fi- 
nal kinetic energy of the fragments with collision velocity ex- 
hibit a more quantitative description. The ability to capture these 
transitions justifies the design of the four-population model and 
demonstrates its descriptive power 

Despite its narrow parameter range, we can draw the follow- 
ing conclusions from the velocity study of Sec. |5] For the transi- 
tion between bouncing and sticking, we find a threshold velocity 
of ~ 1 ms"' and that the transition from sticking to fragmenta- 
tion lies between 11. 5 and 12.5 ms"'. Since our initial filling fac- 
tor (p - 0.35 is close to the critical filling f actor 0.4, wh i ch sep - 
arates porous from compact aggregates in IGiittler et al.l (l20IOl) . 
our findings cannot be directly compared to their results. For 
collisions between equally sized highly porous and very com- 
pact aggregates, they found a direct transition from bouncing to 
fragmentation at ~ 1 ms"'. However, for collisions between a 
compact target and a porous projectile of equal size they find 
a bouncing-sticking transition at 1 ms"' and a transition from 
sticking to no mass gain at 9.4 ms"'. These thresholds resem- 
ble our findings for intermediate porosity very well. Our re- 
sults for fragmenting collisions indicate that the power-law in- 
dex of the fragment mass distribution is velocity dependent. This 
is supported by the collection of laboratory fragmentation data 
( Mathis et al. 1971, iDavis & Rvar]ll990L iBlum & Miin"ci]ll99l 
Giittler et al.,.2010i) . Furthermore, it is physically reasonable that 
in more violent collisions the fraction of smaller fragments in- 
creases. This suggests a velocity-dependent power-law index. 

We note that the results presented in Sec. |5]are valid for col- 
lisions between homogeneous aggregates of intermediate poros- 
ity. Increasing inhomogeneity might affect the presented thresh- 
old velocities. Furthermore, rotation of the largest and second 
largest fragments might cause them to fall apart beyond the sim- 
ulated time, an outcome that eventually affects their final size. 
The results presented here, in particular for fragmenting colli- 
sions, are based on one simulation for each collision velocity. 
Although we do not expect a large variation in outcome because 
of the symmetry given by a head-on collision of two spheres, 
a profound statistical investigation has to be carried out. It re- 
mains to be investigated whether the deviation of the power-law 
mass distribution from a power-law for low fragment masses is a 
numerical or physical effect. Despite these drawbacks, we have 
been able to demonstrate the applicability and functionality of 
the four-population model by means of our simulation results. 

Despite increased experimental efforts, only small spots of 
the required parameter space of pre-planetesimal collisions are 
actually covered by empirical data. Vast regions of these maps 
remain terra incognita and collision data for pre-planetesimal 
sizes larger than centimetre are missing. 

With both a code calibrated for the simulation of pre- 
planetesimals (IGeretshauser et all2010l) and the four-population 
model as an adequate mapping model, we have established a ba- 
sis for profoundly investigating all aspects of pre-planetesimal 
collisions and transferring acquired data to global dust coagu- 
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Table A.l. List of symbols 



symbol 


explanation 


b 


impact parameter 




rotational energy of the target 


P 


rotational energy of the projectile 




second irreducible invariant 


K{<l>) 


bulk modulus 


Ko 


bulk modulus for the RBD dust sample 


m 


mass 


^cum 


cumulated mass 


nif 


fragment mass 


mp 


projectile mass 




mass of the total power-law fragment population 




mass of the total sub-resolution population 


Ml 


target mass 


mi 


mass of the largest fragment 


m2 


mass of the second largest fragment 


P 


hydrostatic pressure 


Pm 


mean pressure of S 


S aft 


deviatoric stress tensor 


T{<P) 


tensile strength 


Va 


velocity vector 


vo 


relative collision velocity 




position vector 




shear strength 


A 


slope of S 


K 


index of the power-law mass distribution 




shear modulus 




most massive member of the power-law population 


P 


density 




stress tensor 


filling factor 




threshold filling factor from elastic to plastic compression 


rc 


threshold filling factor from elastic to plastic tension 




maximum filling factor of S 




minimum filling factor of S 


ipf 


filling factor of the projectile 


4'mu 


filling factor of an uncompressed RDB dust sample 


'Pi 


filling factor of the target 


<t>i 


filling factor of the largest fragment 


02 


filling factor of the second largest fragment 


S(0) 


compressive strength 



lation models. In future studies, these tools can be utilised to 
generate a catalogue of pre-planetesimal collisions. The four- 
population model will be applied to investigate how the collision 
behaviour depends on important parameters such as aggregate 
porosity and inhomogeneity, mass ratio of the collision partners, 
impact velocity, impact parameter, and rotation. Furthermore, 
we will assess the statistics of our fragmentation results. 
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